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Abstract: Multi-objective non-linear programs occur in various field of engineering application. One of the 
applications of such program is structural design problem. In this paper, we consider a generalized form of a 
multi-objective structural design problem. Triangular norm based fuzzy programming technique is used to solve 
these problem. The test problem includes a three-bar planar truss subjected to a single load condition. The model 
is illustrated with numerical examples. 
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I. INTRODUCTION 

Optimization is the process of minimizing or maximizing an objective function (e.g. cost, weight) of a 
structural system which has been frequently employed as the evaluation criterion in structural engineering 
applications. But in the practical optimization problems, usually more than one objective are required to be 
optimized, such as minimum mass or cost, maximum stiffness, minimum displacement at specific structural 
points, maximum natural frequency of free vibration, and maximum structural strain energy. This makes it 
necessary to formulate a multi -objective optimization problem. The first note on multi -objective optimization 
was given by Pareto; since then the determination of the compromise set of a multi -objective problem is called 
Pareto optimization. That is why the application of different optimization technique [11,16-19] to structural 
problems has attracted the interest of many researchers. 

In conventional mathematical programming, the coefficient or parameters of mathematical models are 
assumed to be deterministic and fixed. But, there are many situations where they may not be exactly known i.e., 
they may be somewhat uncertain in nature. Thus the decision making methods under uncertainty are needed. 
The fuzzy programming has been proposed from this point view. In decision making process, first Zadeh[2] 
introduced fuzzy set theory. Tanaka et al.[20]applied the concept of fuzzy sets to decision making problems by 
considering the objectives as fuzzy goals. Later on Bellman and Zadeh [3] used the fuzzy set theory to the 
decision making problem. Zimmermann [4] proposed a fuzzy multi -criteria decision making set, defined as the 
intersection of all fuzzy goals and their constraints. 

In practical, the problem of structural design may be formed as a typical non-linear programming 
problem with non-linear objective functions and constraints functions in fuzzy environment. Some researchers 
applied the fuzzy set theory to Structural model. For example Wang et al. [1] first applied a -cut method to 
structural designs where the non-linear problems were solved with various design levels a , and then a 
sequence of solutions are obtained by setting different level-cut value of a . Rao [8] applied the same a -cut 
method to design a four-bar mechanism for function generating problem .Structural optimization with fuzzy 
parameters was developed by Yeh et.al [7]. In 1989, Xu [6] used two-phase method for fuzzy optimization of 
structures. In 2004, Shih et.al [9] used level-cut approach of the first and second kind for structural design 
optimization problems with fuzzy resources .Shih et.al [10] develop an alternative a -level-cuts methods for 
optimum structural design with fuzzy resources in 2003. Dey et.at [5] optimize structural model in fuzzy 
environment. 

Alsina et.al. [13] introduced the t-norm into fuzzy set theory and suggested that the t-norms be used for 
the intersection of fuzzy sets. Different types of t-norms theory and their fuzzy inference methods were 
introduced by Gupta et.al. [14] .The extension of fuzzy implication operators and generalized fuzzy methods of 
cases were discussed by Ruan et.al. [15]. 

In this paper we propose an approach to solve multi-objective structural model using t-norms based 
fuzzy optimization programming technique. In this structural model formulation, the objective functions are the 
weight of the truss and the deflection of loaded joint; the design variables are the cross-sections of 
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the truss members; the constraints are the stresses in members. The test problem includes a three -bar planar truss 
subjected to a single load condition. This approximation approach is used to solve this multi -objective structural 
optimization model. 

The remainder of this paper is organized in the following way. In section II, we discuss about structural 
optimization model. We discuss about mathematics Prerequisites and aggregation operator in section III and IV 
respectively. In section V, we discuss fuzzy optimization technique to solve multi-objective non-linear 
programming problem. In section VI, we discuss Pareto optimality test. In section VII, we solve multi-objective 
structural model using t-norms based fuzzy optimization. In section VIII, numerical solution of structural model 
of three bar truss. Finally we draw conclusions in section IX. 

IT. MULTI-OBJECTIVE STRUCTURAL MODEL 

In the design of optimal structure i.e. lightest weight of the structure and minimum deflection of loaded 
joint that satisfies all stress constraints in members of the structure. To bar truss structure system the basic 
parameters (including the elastic modulus, material density, the maximum allowable stress, etc.) are known and 
the optimization's target is that identify the optimal bar truss cross-section area so that the structure is of the 
smallest total weight, the minimum nodes displacement, in a given load conditions. 
The multi-objective Structural model can be expressed as: 

Minimize WT(A) 

minimize S(A) 

subject to a{A) < [<r 0 ] 
A min < A < A m . dx 

where A = [Aj,A 2 , ,A n 1 ^ are design variables for the cross section, n is the group number of design variables 

for the cross section bar, WT = ^p^Lj is the total weight of the structure, S(A) is the deflection of loaded 
joint Lj , Aj and p i were the bar length, cross section area, and density of the i' h group bars respectively. cr(A) is 
the stress constraint and[cr 0 ]is maximum allowable stress of the group bars under various conditions, A min and 
are the minimum and maximum cross section area respectively. 

III. PREREQUISITE MATHEMATICS 

111(a). Fuzzy Set: 

Let X is a set (space), with a generic element of X denoted by x , that is X(x) .Then a Fuzzy set (FS) is 
defined as A = {(x,fi A (x)) : x e x) 

where ju^ : X — > [0,1] is the membership function of FS £4 . ju^ (x) is the degree of membership of the element 
x to the set 54 . 

111(b). a -Level Set or a -cut of a Fuzzy Set: 

The a -level set of the fuzzy set 14 of X is a crisp set A a that contains all the elements of X that have 
membership values greater than or equal to a i.e. A = : (x) > a, x e X, a e [0, 1]J . 
III(c). Convex fuzzy set: 

A fuzzy set 14 of the universe of discourse X is convex if and only if for all x { ,x 2 in 
X, / / J (Ax 1 +(l-Ax 2 )>min(// a (xi),// a (x 2 ))) whenO<A<l. 

IV. AGGREGATION OPERATOR 

When the rules in the decision support system contain more than one antecedent, the degrees of 
strength of antecedents need to be combined to determine the overall strength of the rule consequent. In the 
language of fuzzy sets, the membership values of the linguistic variables in the rule antecedents have to be 

combined using an aggregation operator. Formally, a general aggregation is a real function T : [0,1]" — >[0,l] , 

non decreasing in all arguments, with the properties 7(0) = 0 and T(l) = 1 . 

General aggregation operators display the whole range of behavior, disjunctive, conjunctive, averaging, mixed, 
commutative, mutually reinforcing or otherwise and correspond to vague and loosely defined "and" and "or" 
connectives etc. Triangular norms and conforms and averaging operators are well known examples of the 
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aggregation operators. Different class of aggregation operators display substantially different behavior, it is not 
logical to use any particular class to provide generic representation of aggregation. Therefore, we will use 
general aggregation operators to model aggregation of rule antecedents in decision support systems. They will 
provide the highest degree of adaptability and excellent empirical fit. However, if there are strong reasons to 
restrict the selection to a particular family of operators, we will impose the relevant constraints. 
Consider general aggregation operator. The function can have a simple algebraic form, such as 
T(x 1 ,x 2 ,....,x n ) = wm{x 1 ,x 2 , x„} 



The degree of importance of rule antecedents (vector a ) can be easily incorporated into aggregation operators 
in variety of ways. For example 

T(x l ,a l ;x 2 ,a 2 ; ;x n ,a n ) = min{x 1 ,a 1 }xmin{x2,a 2 } x xmm {-*n> a n} 

T(x l ,a 1 ;x 2 ,a 2 ; ;x„,a„) = min{x 1 a 1 +x 2 a 2 + + x n a n ,l) 

In this article, decision making method used by the weighted bounded sum operator (member of Yager family of 
triangular conorms). 

V. FUZZY PROGRAMMING TECHNIQUE TO SOLVE MULTI-OBJECTIVE 

NON-LINEAR PROGRAMMING PROBLEM (MONLP) 

A Multi-Objective Non-Linear Programming (MONPL) or Vector Minimization problem (VMP) may 
be taken in the following form: 

Minimize f(x) = [f x (x), f 2 (x), f k (x)f (2) 

subjectto xeX={xefl" : gj (x)<or =or>bj for j =1,2,3,. ...,m\ and /, < x, < u t (i = 1,2,3, n) 

Zimmermann (1978) showed that fuzzy programming technique can be used to solve the multi -objective 
programming problem. 

To solve MONLP problem, following steps are used: 
Step 1: Solve the MONLP (2) as a single objective non-linear programming problem using only one objective at 
a time and ignoring the others, these solutions are known as ideal solution. 

Step 2: From the result of step 1, determine the corresponding values for every objective at each solution 
derived. With the values of all objectives at each ideal solution, pay-off matrix can be formulated as follows: 





A(x) 


/ 2 W 


f k (x) 


x 1 




/ 2 V) ■ 




x 2 


/iV) 


/ 2 V) • 


■■ K(x 2 ) 


x l 


/iV) 


/ 2 V) 


AV) 






J 



Here x x ,x , , x are the ideal solutions of the objectives /i (x) , f 2 (x) , , f k (x) respectively. 

So U r =rmx[f r {x l \f r {x 2 \ ,f r (x k )} and L r =min{/,.(x 1 ),/ r (x 2 ), ,/ r (x*)} 

Where t/ r and L r be upper and lower bounds of the r"' objective function f r (x) for r = 1,2, 3, , k. 

Step 3: Using aspiration level of each objective of the MONLP (3) may be written as follows: 
Find x so as to satisfy 

f r (x)< L r with tolerance P r (=U r -L r ) for r = 1,2,3, , k 

xeX. l.^x^Ui (i = 1,2,3,. ...,«) 
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Here objective functions of (2) are considered as fuzzy constraints. These types of fuzzy constraints can be 
quantified by eliciting a corresponding membership function: 



(3) 



[1 if f r (x)<L r 

where L r = L,.+ e,.and 0 <e r < U r -L r , for r = 1,2,3,..., A: 

Mr(fM)) 




Figure 1: Membership function for objective functions /,.(.v) 
Having elicited the membership functions as in (3) ju r (/,. (x)) for r = 1,2,3,..., & a general aggregation function 

jLfy (x) = F( Ml (/i (x)),m 2 if 2 (*))- >Mk (fk ))) is introduced. 

So a fuzzy multi -objective decision making problem can be defined as 

Maximize 
subject to x e X. 

Fuzzy decision making method used by the (weighted) bounded sum operator (member of Yager family of 
triangular conforms) the problem (4) is reduced to 

Maximize ju$ (x;W) = ^JV r M r (f r (x)) 

subject to xe X 

(5) 

0<ju r (f r (x))<l for r = \,2,~,k 

where W r > 0 for all r = \,2,...,k, ^JV r =1 

Step 4: Solve (5) to get optimal solution. 

Some basic definitions are introduced below. 
V(a). Complete Optimal Solution 

x is said to be a complete optimal solution to the MONLP (3) if and only if there exists xeX such 
that/ r (x* j < f r (x) for r = \,2,...,k and for allxeX .However, when the objective functions of the MONLP 

conflict with each other, a complete optimal solution does not always exist and hence the Pareto Optimality 
Concept arises and it is defined as follows. 

V(b). Pareto Optimal Solution 
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x is said to be a Pareto optimal solution to the MONLP (3) if and only if there does not exist 
another x e X such that f r [x ) < f r (x) for all r = 1,2,..., k and f j (x)*f j {x*^ for at least one j, 
ye {1,2,...,*}. 

VI. PARETO OPTIMALITY TEST 

A numerical test of Pareto optimality for X can be performed by solving the following problem: 



Maximize R = y e, 



(6) 



subject to f r (x) + e r = f r [x ) ,r = 1, 2, ...,k 
xeX; s r >0. 

The optimal solution of (6), say x** and f r [x' ^ are called strong Pareto optimal solution provided V is very 
small otherwise it is called weak Pareto solution. 

VII. FUZZY PROGRAMMING TECHNIQUE IN MULTI-OBJECTIVE 

STRUCTURAL MODEL 

To solve the above MOSOP (1), step 1 of V is used. After that according to step 2 pay-off matrix 
formulated as follows: 

WT(A) S(A) 



WT\A 1 *) 5{A 1 *) 
WT(A 2 *) 8"(A-) 



After that according to step 2, the bounds of objective are U l , L, for weight function WT(A) 

(where L, < WT( A) <U l ) and the bounds of objective are U 2 ,L 2 for deflection function 8(A) 

(where L, < S(A) < U 2 ) are identified. 

Above MOSOPP reduces to a FMOSOPP as follows; 

Find A 

Such that 

WT(A)< Lj with maximum allowable tolerance P { |= U l -Ly j 
<5( A)< L 2 with maximum allowable tolerance P 2 U 2 — L^j 
a(A)<[a Q ] 

Aran - ^ - Anax 

Here for simplicity linear membership functions ju^^WT^A)) and fi s (S(A)) for the objective functions 
Wr(A) and S(A) respectively are defined as follows: 
[l ifWT(A)<L l 
( U X -WT(A)\ 



MwT {WT(A)) = 



U.-L, 



M S ( A )) = 



z/L, <WT(A)<U 1 

if WT(A)>U l 
ifS(A)<I^ 

if L 2 <S(A)<U 2 

if 5(A)>U 2 
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According to step-3, having elicited the above membership functions crisp non-linear programming problem is 
formulated as follows 

Maximize {w i fi wr (wr(A))+W 2 ju s (S(A))) ^ 
subject to 

0< ju wt (WT(A))<1, 0<ju s (S(A))<l, 
er(A) < [a], 
A ran < A< A max , 
W l >0,W 2 >0, W l +W 2 =1; 
The problem (7) can be written as 

■Ju.-mW^fUz-siA)]) 



sub/eel 1 



( U l -WT(Ay 



<1, 0< 



If-SiA)) 



a(A) < [a], 
A^ < A< A max , 
W, >0,W 2 >0, Wi+W 2 =1: 



U 2 -L 2 



VIII. NUMERICAL SOLUTION OF A MULTI-OBJECTIVE STRUCTURAL 

OPTIMIZATION MODEL OF A THREE-BAR TRUSS 

A well-known three bar [12] planar truss structure is considered. The design objective is to minimize 
weight of the structural WT(A X ,A 2 ) and minimize the deflection S(A l ,A 2 ) along x and y axes at loading 
point of a statistically loaded three -bar planar truss subjected to stress (cr) constraints on each of the truss 
members. 




Figure 2. Design of the three-bar planar truss 

The multi-objective optimization problem can be stated as follows: 
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Minimize WT(A l ,A 2 ) = pL(2A l +A 2 ) 

. PL(2A l +A 2 ) 
minimize 8 X [ A l , A 2 ) = 

minimize S y (Ay, Az) 



£(2 A, 2 + 2A 1 A 2 ) 

PLA 2 
e(2A}+2A,A 7 \ 



x P(2A,+A 2 ) r T -, 

subject to o-! ( 4 , A 2 ) = ~^f^ ~ [°i ] (9) 

^^'^^^ 

ct 3 (A 1 ,A 2 ) = ^a—^rofl 

~ V 7 2^2+2^2 L J 

A mn <A<A max ;' = U 

Where P = applied load; p =material density, L ^Length, E=Young's modulus, A l = cross section of bar-1 and 
bar-3, A 2 =cross section of bar-2. 8 r and 8 are the deflection of loaded joint along x and y axes 
respectively, [of] an d \_ a 2~\ me maximum allowable tensile stress for bar 1 and bar 2 respectively, [erf J is the 
maximum allowable compressive stress for bar 3. 

The input data for MOSOP (9) is given as follows 



Appli 

ed 
load 

P 
(KN) 


Volume 
density p 


Length 
L 

(«) 


Maximum 
allowable 
tensile 
stress 

K] 

(KN/m 2 ) 


Maximum 
allowable 
tensile 
stress 

(KN/m 2 ) 


Maximum 
allowable 
compressive 
stress 

[•?] 

(KN/m 2 ) 


Young's 
modulus E 

(KN/m 2 ) 


A™ n and 
A™" of 
cross section 
of bars 

(W) 


20 


100 


1 


20 


10 


20 


2xl0 8 


A™" =0.1 
A, max = 5 
AT" =0.1 

A™ x = 5 



Solution: According to step 2 pay off matrix is formulated as follows: 

8 x (A l ,A 2 ) 



f 2.187673 
15 
10.1 

V. 



<?,(A„A 2 ) 
5.857864 



3 

3.960784 



J 



Here U m =15 , 1^=2.187673, £^=3 E/^=20, £^=3,^=3.5, [/^ =5.857864, L s =0.039216, 
L s =0.15; Here linear membership for the objective functions WT(A i ,A 2 ) , ^(A^A^ and 8 (Ay, A 2 )are 
defined as follows: 
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[1 ifWT(A l ,A 2 )<3 
Mm (WT(A, , A 2 )) = 15 ~ WT (Ai,A 2 ) g. 3 ^ , A 2 ) < 15 

[o if WT(A,,A 2 )>15 




3 15 
Figure 3. Membership and non-membership function for WT(A l ,A 2 ) 

fl if S X {\,A 2 )<3.5 



if S X (A U A 2 )> 20 



//(5 X (A,A 2 )) 




Figure 4. Membership and non-membership function for S^^Aj) 

fl i/ ( J V (A 1 ,A 2 )<0.14 

5.857864 -<?, (A, ,A 2 ) 

H S (SJA l ,A 2 )) = \ } - — i/0. 14 <^ V (A!,A 2 )< 5.857864 

™' y 1 2 5.717864 - 1 1 2 

0 ;/ 5 (A 1 ,A 2 )> 5.857864 
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ju(S y (A,Aj) 




0.14 5.857864 
Figure 5. Membership and non-membership function for 5 (4,4) 
Now Fuzzy decision making method used by weighted bounded sum operator (member of Yager family of 
triangular conorms) , 

Maximize F = W l ju(WT(A l ,A 2 )) + W 2 ju(SJA l ,A 2 )) + W } ju(S Y (A l ,A 2 )) (10) 
subject to 

0<fi(WT(A l ,A 2 ))<l; 0<ju(S x (A i ,A 2 ))<l 
0<p(S y (A 1 ,A 2 ))<l 
15-WT(A 1 ,A 1 ) 



J u(WT{A l ,A 2 )) = 



/,(<? t (A p A 2 )) 



12 

20-^(4,4) 



16.5 

5.857864-^(4,4) 
5.717864 

20(24 + 4) ^ 
24^244 " 

20 



%/2(4+4) 

20(A 2 ) 

24^244 
o.i<4,4 <5, 

W,+W 2 +W 3 =l, 



: 20: 
C10; 
C20; 



The solution obtained from Eq. (10) is given in table 2. 

Table 2: Optimal results of MOSOP (9) 



w l 


W 2 


^3 


A* 


4* 


WT*(4,4) 


<(4,4) 










xl0~ 4 m 2 


xl0~ 4 m 2 


x\0 2 KN 


xl0" 7 m 


xl0" 7 m 


1/3 


1/3 


1/3 


2.661308 


0.1029932 


5.425610 


7.375100 


0.14 


0.6 


0.2 


0.2 


1.645225 


0.1 


3.390449 


11.80812 


0.3482759 


0.2 


0.6 


0.2 


4.542762 


0.3 


9.394060 


4.262608 


0.14 


0.2 


0.2 


0.6 


2.661309 


0.1029933 


5.425611 


7.375099 


0.14 



In table 3, the value of R is quite small and hence the optimal results in table 2 are strong Pareto-optimal 
solution and can be accepted. 



International organization of Scientific Research 



42 I P a g e 



Fuzzy Optimization Technique for Pareto Optimal Solution of Structural Models with Stress constraints 



Table 3: Pa reto optimal results of MOSOP (9) 



R 


A" 


A" 


WT**(A l ,A 2 ) 


<C(AA) 








xl0~ 4 m 2 


xlO 2 KN 


xl0~ 7 m 




0.5277774 xlO" 6 


2.661308 


0.1029933 


5.425609 


7.375100 


0.14 



IX. CONCLUSIONS 

The present paper proposes a solution procedure for structural model. The results of this study may 
lead to the development of effective t-norm optimization method. This solution procedure may be used for 
solving other model of nonlinear programming problem in different field. 
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